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; Abstract 

■ To obtain a well defined path integral one often employs discretizations. In the case of 
gravity and reparametrization invariant systems, the latter of which we consider here as a toy 

ly^ ■ example, discretizations generically break diffeomorphism and reparametrization symmetry, 

] respectively. This has severe implications, as these symmetries determine the dynamics of 

the corresponding system. 

^ ' Indeed we will show that a discretized path integral with reparametrization invariance 

<sJ '< is necessarily also discretization independent and therefore uniquely determined by the cor- 

^ ' responding continuum quantum mechanical propagator. We use this insight to develop an 

iterative method for constructing such a discretized path integral, akin to a Wilsonian RG 
flow. This allows us to address the problem of discretization ambiguities and of an anomaly- 
free path integral measure for such systems. The latter is needed to obtain a path integral, 
^ ' that can act as a projector onto the physical states, satisfying the quantum constraints. We 

. will comment on implications for discrete quantum gravity models, such as spin foams. 

■ 1 Introduction 

Discretizations have become a popular tool in classical and quantum physics, on the one hand to 
allow for a numerical treatment, on the other hand to regularize for instance the path integral. 
Also a quite popular expectation is that fundamental physics is based on discrete structures 
rather than continuum space time, a conjecture followed in a number of quantum gravity ap- 
proaches, for instance [Hll]. 
5^ I However, in discretizing a continuum theory one has to face several issues. One issue, which 

is particularly relevant for discrete quantum gravity approaches, is that the symmetries of the 
continuum theory might be broken by the discretization. Another issue is that discretizations 
are typically never unique, many different discrete models might lead to the same continuum 
physics This is not a problem if discretization is just viewed as a tool for obtaining continuum 
physics, but has to be addressed if the discrete theory is claimed to be 'fundamental'. 

These problems appear in one form or other in many approaches to quantum gravity, where 
discretizations break diffeomorphism symmetry [3] . This symmetry is deeply entangled with the 
dynamics of the theory and hence its breaking has particularly severe repercussions [Ij. 



X 



^Although one of the main problems in many discrete quantum gravity approaches is actually to extract large 
scale physics. 
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In the canonical framework diffeomorphism symmetry leads to the Hamiltonian and diffeo- 
morphism constraints, satisfying the so-called Dirac algebra, a canonical version of the group 
property of diffeomorphisms. These constraints are central to the (quantum and classical) dy- 
namics. A discretization, violating diffeomorphism symmetry leads however to a violation of the 
Dirac algebra and inconsistencies in the dynamics [5l[6]. Some of these issues can be addressed 
in the consistent and uniform discretization [7] and master constraint [8J approach, however 
there the symmetries can only be regained in the continuum limit (if at all) [9J. There is an 
anomaly-free quantization of the Hamiltonian constraints [T^ in Loop Quantum Gravity, how- 
ever it features many ambiguities whose significance is not fully understood yet [llj and the 
status of the corresponding Dirac algebra is not fully satisfactory \10\ [12] . Although one might 
suspect that in a covariant approach the situation regarding diffeomorphism symmetry is much 
better than in the canonical ones, this is not the case - indeed the problems are closely con- 
nected [21 m [T3l I14j . So far, we are lacking a path integral in which diffeomorphism symmetry 
is realized on the discrete level. Related to this issue is the lack of 'discretization independent' 
(four-dimensional gravity) models, as there are arguments [T5] that discretization independence 
should be the equivalent to the diffeomorphism invariance of the continuum. Indeed, we will 
present in section [2] an argument involving the dynamics of the theory, showing this equivalence 
for the toy model of a one-dimensional reparametrization invariant system. 

A proper (non-compact) gauge symmetry in a path integral leads to divergencies, as one 
integrates an invariant amplitude over gauge orbits. If the symmetries are broken, these diver- 
gencies do not necessarily appear on the discrete level, but should reappear in the continuum 
limit. However, these divergencies have to be distinguished from proper UV divergencies which 
might also appear in a continuum limit. Additionally in a perturbative approach to path inte- 
grals broken symmetries might lead to inconsistencies |14] . 

In addition to the amplitudes for the path integral one has to define an integration measure. 
Usually one would like this integration measure to be anomaly-free, i.e. to respect the gauge 
symmetries of the action. Although there is no discrete amplitude respecting diffeomorphism 
invariance yet, there is some work [16 [ I17 [ [TBj deriving invariance conditions for the measure in 
spin foam models and in this way fixing some of the ambiguities. These however do not involve 
the full dynamics so far. Again we will show in the context of the one-dimensional toy model 
that a symmetry preserving measure can be found. For this one has however to involve the full 
dynamics of the theory. 

There is a wide range of ideas to tackle these issues, which we cannot all discuss here. We 
will rather illustrate one approach, which is to obtain so-called perfect discretizations via a 
Wilsonian renormalization procedure [19[ \20\ \2T\ [22] , within the toy model of a one-dimensional 
reparametrization invariant system. The basic idea is to 'pull back' continuum physics onto the 
lattice. Alternatively one can start with a discretization, subject it to a renormalization group 
transformation, which maps the physics of a fine discretization to a coarse grained one. Iterating 
this procedure one should be able to find a fixed point which gives the continuum physics exactly 
mirrored onto the discretization. 

Although the model we will consider in this work is extremely simple compared to gravity, 
it displays many features that we believe will also be relevant for gravity. We will show that 
even for this simple system, a perfect discretization cannot be simply guessed, but has to be 
determined by solving the dynamics of the system. In particular we will show that a full 
implementation of reparametrization invariance - the equivalent to diffeomorphism invariance 
in general relativity - into the discrete path integral implies 'discretization invariance', i.e. the 
result of the path integral is independent of the number of subdivisions, including the possibility 
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to have none at all. We expect that a similar feature will hold for gravity. This 'discretization 
invariance' actually means that the basic discrete amplitude (and measure) provides already the 
full quantum propagator of the system, also on large scale. Hence any free parameter in such 
a discretization will influence macroscopic physics. In this sense, requiring reparametrization 
invariance for the discretization, will resolve all discretization ambiguities. We can expect that 
likewise many ambiguities for discrete quantum gravity are fixed by requiring diffeomorphism 
invariance. The crucial question is, whether we are left with finitely or infinitely many (relevant) 
parameters [23] . 

As will be explained in section[2]for the one-dimensional reparametrization invariant systems, 
a perfect discretization of the path integral is provided by the (continuum) quantum propagator, 
i.e. the evaluation of the continuum path integral. This is however not very helpful in the many 
situations where such a continuum quantum propagator or path integral is not available. In this 
case one typically rather has to start with the discrete system and develop some approximation 
methods. We therefore develop an iterative method to evaluate the path integral, basically 
following a Wilsonian renormalization approach. This approach will automatically address the 
issue of discretization ambiguities, as we do not consider a specific discretization but a parameter 
space describing a certain class of discretizations, in which one determines the flow induced by 
a renormalization group transformation. 

A technical difficulty which we address here, is to deal with the broken gauge symmetries 
in the path integral. This can be done by expanding around a special solution, for which the 
symmetries are exact. (Typically the zero solution = for one-dimensional systems. For 
many discretizations of gravity flat space is such a solution.) The action (and measure) can 
then be improved order by order such that the gauge symmetries are exact to the order where 
one has to evaluate the path integral. Hence one can either gauge flx or cleanly separate the 
(infinite) factor arising from the integration over the gauge orbits. 

To start with we will explain the basic ideas of the approach in the next section [2l In 
particular we will show that for one-dimensional systems reparametrization invariance and dis- 
cretization independence for the discrete path integral are equivalent. In the following section [3] 
we will consider the parametrized, discretized harmonic oscillator, and demonstrate the concepts 
from section [2j Since in this case the time variable is discrete as well, this is not a linear system 
anymore, and we consider different linearizations in order to compute the reparametrization- 
invariant propagator. In section [5] we will treat the (quartic) anharmonic oscillator with the 
same methods, computing the perfect propagator to first order in the interaction parameter 
A. Finally, we will summarize and discuss our findings in section [5l considering in particular 
the implications for discrete gravity and spin foam models. The appendices El [Bl O contain 
more technical details on the uniqueness of the solutions to the recursion relations describing 
the renormalization flow, and techniques to solve such recursion relations. 

2 Reparametrization invariance in the discrete and discretiza- 
tion independence 

Consider a one-dimensional mechanical system defined by a Lagrangian L{q, q). The variation of 
the corresponding action S = j L{q, q) with respect to q{t) will determine the solutions. Such a 
system can be made reparametrization invariant by adding the time parameter t as a dynamical 
variable to the configuration variable q. Time evolution is then with respect to an auxilary 
parameter s. One can define a new action (here ' denotes the derivative with respect to the 
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auxiliary parameter s) 

S{qis),t{s)) = I dsLl^q,^^' (2.1) 

which is invariant under reparametrizations q{s),t(s) — )■ q{f{s)),t{f{s)) of the trajectories. This 
new action has to be varied with respect to q{s) and t{s). The equation of motion for t{s) wih 
however be automatically satisfied, if the equation of motion for q{s) is. Hence our system is 
underdetermined as we have two variables but only one independent equation of motion. Indeed, 
the solutions (?s(s), ts{s) are not uniquely determined by the boundary conditions, as given one 
such solution one can find a family of (physically equivalent) solutions Q's(/(s)), is(/(s)) by 
reparametrizing the auxiliary evolution parameter s — t- f{s). Performing a (singular) Legendre 
transformation one will find instead of a proper Hamiltonian a Hamiltonian constraint C = pt+H 
(where H is the Hamiltonian of the original system), which at the same time is the generator of 
the gauge transformations, i.e. it generates evolution in the auxiliary parameter s. This is very 
similar to general relativity (where the role of the auxiliary parameter is taken over by space 
time coordinates). 

This reparametrization symmetry is typically broken if we discretize the system. To this 
end we replace q{s),t{s) with s in some finite interval, by some finite set of variables qn,tn with 
n = 0, . . . , N . One method of discretization is to replace derivatives by difference quotients and 
to choose some discretization for the potential term V, for instance 

Af-l Af-1 

S{qn:tn) ■= := ^ 

n=0 n=0 

(Here and in the following we will consider Wick rotated actions in order to make the path 
integrals convergent. This just changes the minus in front of the potential F to a plus sign.) 
The factor {tn+i — tn) arises from the integration measure / dst' in (j2.ip . 

This action (|2.2p does in general (an exception being the free particle V = 0) not feature 
anymore any gauge symmetries. As a result, the equations of motion arising by varying with 
respect to t„, n = 1, . . . , iV — 1 in fact uniquely fix both qn and tn- Hence both q and t become 
propagating (or physical) degrees of freedom in the discrete theory. Only in the continuum limit 
does reparametrization invariance arise, so that t becomes a gauge variable again. The broken 
symmetries lead to so-called pseudo constraints in the canonical formalism [3l \7\ [H] instead of 
proper constraints. These are equations of motion, i.e. equations between the canonical data of 
two consecutive time steps n and n + 1, which do however only weakly depend on the data at 
n + 1. (Proper constraints are equations of motion which do involve only data of one time step 
n.) 

This feature arises also for discretizations of gravity, such as Regge Calculus [3]. In this case, 
the discretization (i.e. the triangulation of space-time) breaks reparametrization-invariance (i.e. 
four-dimensional diffeomorphism invariance) as well. Furthermore, in gravity theories singling 
out an equivalent for t, i.e. identifying the (pseudo) gauge degrees of freedom and separating 
them from the truly physical ones, becomes hideously complicated. 

The breaking of gauge symmetries is a result of the choice of discretization, however. In par- 
ticular, there are discrete actions in ID which exhibit a discrete remnant of the reparametrization 
invariance, which lead to the correct amount of physical degrees of freedom. Such perfect ac- 
tions can be constructed by a refinement process. Starting from a discrete action S{qn,tn), one 
improves it by refining the discretization, solving the equations of motion for the refined degrees 
of freedom, and evaluating the refined action on that solution. By iterating this process, or by 



Qn+i — Qr, 



n+l 



tr, 



+ l^{V{qr:) + V{qn+i)) 



{tn+l-tn) (2.2) 
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directly considering the limit of sending the refined discretization to the continuum, one will 
find the perfect action, which can be shown [3] to be given by Hamilton's principal function for 
the system, i.e. 



Sperfiln, t-n) = SuPFil 



(2.3) 



n 



where SHPFiQi,ti,qf,tf) is Hamilton's principal function (for the continuum system), i.e. the 
action evaluated on the solution q{s), t{s) with boundary conditions {qi,ti,qf,tf). The perfect 
action does display a gauge symmetry for every (inner) discretization point n, which we will call 
vertex translation symmetry, as it is of the form (for finite gauge parameter A) 



i.e. it translates the discretization point in time. Hence all the tn can be seen as gauge param- 
eters. Furthermore the action is 'discretization independent', i.e. Hamilton's principal function 
SHPFiQOjto, QNitN) (being the classical equivalent to the partition function or path integral, 
which encodes all the dynamical information) computed from Sperf{qn,tn) does not depend on 
the number of discretization points used, and coincides with the one computed from the contin- 
uum action Shpf{qo, to, qN,tN) = Sperf{qo,to, qN,tN)- In particular we can choose N = 1. Note 
that to construct this action one needs to basically solve the dynamics of the system. This can 
be seen as a huge disadvantage of the method. However Hamilton's principal function provides 
actually the only discretization featuring vertex translation symmetry [24J. Hence addressing 
the dynamics of the system cannot be avoided if one wants to implement the symmetries of the 
system (or even just ensure that these symmetries appear in the continuum limit). 

Instead of solving the dynamics at once, one can consider various approximations and follow 
an iterative approach of constructing better and better actions. We will follow this idea here 
and see that it is closely related to a Wilsonian renormalization group approach. The perfect 
actions then arise as fixed points of iterative (renormalization group) transformations. Note 
that this also resolves (at least in one dimension) all the discretization ambiguities, as Sp^rf is 
the unique action, displaying (the discrete remnant) of discretization invariance. 

It has been shown that this procedure also works for 3D Regge Calculus with nonzero 
cosmological constant For first steps to dealing with AD gravity in a perturbative set-up, 
see [22]. 

This addresses the classical theory. An open issue is, whether a similar approach will work 
for the quantum theory, here in a path integral approach. To evaluate path integrals analytically 
we have to follow a perturbative approach. For broken symmetries such a perturbative ansatz 
may turn out to be inconsistent however [14] . The problem are exceptional solutions (such as 
qn = 0), which are typically the only ones displaying symmetry under vertex translations (i.e. 
the tn can be chosen arbitrarily, the qn remain zero). Perturbing around such solutions one will 
not find a quadratic term for the time variables, however these time variables will appear in the 
higher order potential terms. This hinders the perturbative evaluation of the path integral. A 
way out, suggested in [lll[22| and also followed up here, is to improve the action perturbatively, 
so that gauge invariance can be obtained order by order. This allows an evaluation of the path 
integral to the corresponding order, either by gauge fixing or by changing to gauge invariant 
variables. 

In addition to the action, a path integral requires an integration measure. In the presence 
of gauge symmetries one would require this measure to be invariant under these symmetries. 






5 



otherwise one will obtain anomalies. Indeed only if the measure is invariant, can the path 
integral serve as a projector onto the states satisfying the constraints (arising in a canonical 
quantization) [13]. We will argue here that, similarly to having to solve the classical dynamics 
to obtain the perfect action, one needs to solve the quantum dynamics to obtain the 'perfect 
measurell and with this a perfect discretization of the path integral: 
Consider a discrete path integral with two time steps 

{Qo,to\q2,t2) ■■= Z{qo,to,q2,t2) ■= j dqidti K{qQ,to,qi,ti)K{qi,ti,q2,t2) (2.5) 

where we summarised amplitude and measure for one discretization step into the discrete prop- 
agator K{qn,tn;qn+i,tn+i)- (We term this the propagator as the perfect discretization is given 
by the quantum propagator of the system.) Now assume that vertex translation invariance ()2.4p 
has been fully implemented into ()2.5p . As the gauge symmetry just translates the time variable 
ti we can gauge fix to an arbitrary value ti = t{ and drop the ti integration. (The Fadeev-Popov 
determinant is trivial in this case.) We obtain 

Ziqo,to,q2,t2) ■■= J dqi K{qo,tQ,qi,t{)K{qi,t{,q2,t2) ■ (2.6) 

As t{ is arbitrary we can consider the limit t{ — )• ^2- Now the interpretation of K{qi,t[, 92,^2) is 
to give the amplitude for a particle to propagate from qi to q2 during the time interval (t2 — ^ 1 ) • 
In the limit t{ — > t2 we should therefore have K{qi,t{, 92, ^2) — ^ ^{qi — ^2)- Hence we will obtain 
that the right hand side of ()2.6p is equal to the discrete propagator K{qo,to,q2,t2) and since t{ 
can be chosen arbitrarily we have shown that 

Kiqo,to,q2,t2) = j dqi K{qo,to,qi,ti)K{qi,ti,q2,t2) ■ (2.7) 

Starting from the assumption, that vertex translation invariance has been realized for the path 
integral ()2.5p , we have shown that the discrete propagator K needs to satisfy ()2.7p , which is the 
usual convolution property for the propagator kernel in quantum mechanics. 

This actually proves that a path integral with vertex translation invariance is also discretiza- 
tion invariant, i.e. it does not depend on the number of subdivisions. Having no subdivisions at 
all ~ which gives just the discrete propagator K - should coincide with having infinitely many - 
which gives the propagator in the continuum. Hence the discrete propagator K{qn, tn, g^+ij ^n+i) 
is given by the usual quantum mechanical propagator for the continuum system. 

On the other hand, if we have a discrete propagator satisfying (|2.7p (where we do not integrate 
over t), then the corresponding path integral is invariant under vertex translations, as one can 
integrate out and reinsert every discretization point (g„, t„). This shows that for one-dimensional 
reparametrization invariant systems, finding a discretization which respects this invariance and 
discretization independenc^ are equivalent. We conjecture that this will also hold for discrete 
gravity, i.e. discretized path integrals respecting diffeomorphism invariance (which in its discrete 
form is also expected to be vertex translation invariance, see the discussion in [22]) should be 
also discretization independent. Again, this also means that in order to construct such path 
integrals one has to consider the dynamics of the discrete (quantum) models. 



^The split of the path integral into amplitude and measure is ambiguous and we will not insist on one particular 
splitting here. 

''where here we understand under discretization independence the property (|2.7p . There we do not integrate 
over the time variable, that is we consider already a gauge fixed version of the path integral. Otherwise the 
integral would be divergent, if vertex translation invariance is realized. 
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To summarize the discussion, a perfect discretization of the quantum mechanical path inte- 
gral would be given by the (continuum) propagator. This would however not be very helpful 
for systems in which this propagator is not known, such as quantum gravity. Therefore in this 
article, we want to adopt the procedure of successively improving the action of a classical dis- 
crete system to the quantized this approach might also be helpful for more complicated 
cases. This will also introduce a method to actually solve the path integral for the corresponding 
continuum system iteratively. 

To improve the discrete propagator iteratively we start from the propagator property ()2.7p 
with a 'naively discretized' propagator This can be taken as 

K^'\qn,tn,qn+i,tn+i) = r?f exp(-i5f ) (2.8) 

where S^^^ is the naive discretized action (|2.2p . We will see that the classical part of the iteration 
equations (which is the part without h dependence) leads to the perfect discretization for the 
action. The initial measure factor rjn ^ should be chosen such that K^'^^ — > 5(g„ — q-n-i) for 
(t„_i - tn) 0. 

Refining the discretization and integrating out the intermediate degrees of freedom will 
result in the improved propagator K^-^\ and iterating the procedure will lead to the perfect 
propagator satisfying ()2.7p . From the construction it should be clear that one subdivision step 
will be sufficient, i.e. we define 

^^"+'Hgo,io,'?2,t2) := J dqidhK^^Hqo,to,qi,h)K'^''\quh,q2,t2) . (2.9) 

Here we included the integration over ti . In section [3] we will perturb around the solution 
qn = which is reparametrization invariant also in the naive discretization and introduce a 
method such that the ti integration can be dropped. 



K^''Kqo,to,quti) K(^\qi,ti,q2,t2) 




7f W 



90,^0 92,^2 

K(^+^\qo,to,q2,t2) 

Figure 1: The refinement process of the propagator K^""^ involves subdividing the discretization 
intervals and integrating over the new variables, obtaining a new propagator 

^(n+i)^ This 

process can be iterated, leading in the limit to the perfect propagator, which is in particular 
invariant under refinement of the discretization. 



In order to find the perfect propagator, we will parameterize (some part of) the space of 
functions, allowing for arbitrary couplings in the action and measure factors. Equation (|2.9p 
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will then lead to recursion relations for the coefficients, which are closely related to the renor- 
malization group flow of the system. We will find that, not only does the classical perfect action 
provide a fixed point of the equations derived from ()2.9p . but also the measure factor r/ will be 
determined, providing the quantum corrections to the classical perfect action. 

3 The harmonic oscillator linearized around g„ = 

In the following, we will deal with the (Wick-rotated) parameterized harmonic oscillator, the 
action for which is given by 

S{q{s), t{s)) ds ^'^+ {s)^ . (3.1) 

By naively discretizing the action one replaces continuum functions q{s), t(s) by discretely 
many values qn = q{sn), tn = t{sn), and derivatives by difference quotients. For the discretization 
of the integral measure we use ds = Sn+i — Sn- Hence we arrive at 

n=0 

As we have already indicated, because of the naive discretization, the action (j3.2p does not con- 
tain gauge symmetries, see for instance [3]. That is given (general) boundary data {qo,to;qN,t]y) 
both the qn and the t„ for n = 1, . . . — 1 are determined by the equations of motion. 

The action (|3.2p is not quadratic and not even polynomial in the t„-variables, and therefore 
not considered to be a non-interacting system. The path integral is therefore non-trivial, and 
we consider approximations to it which will render it feasible. The approximation we are going 
to perform here involves quantizing only the perturbations around a background solution, i.e. 
we consider the perturbed variables 

qn = qn + Xn (3.3) 
tn = tn + Tn (3.4) 

where g„, in are to satisfy the equations of motion coming from the discrete action. The path 
integral for this linearized theory, containing the variables x„, can be viewed as a linear 
approximation to the path integral for the whole theory. Of all the solutions to the equations 
of motion derived from (|3.2p . the constant solution, i.e. qn = 0, tn arbitrary, is special in the 
regards that this solution still exhibits some gauge symmetry (i.e. change in the tn), which is 
a remnant of the reparametrization invariance of the continuum theory. In this regards, it is 
analogous to the flat solution et = in AD Regge calculus [3l [26] , which still exhibits the vertex 
displacement symmetry. We will consider linearization around this solution first, before turning 
to linearization around a general solution later. 

Expanding the action for one time-step (neglecting all orders of x„, r„ higher than two), we 
arrive at 

5W(rro,ro,.i,n) = ^^^i-^ + ^{xj + xDn, (3.5) 
Z ioi 4 

=: af\Toi){xl + xj) + af\Toi)xoXi 



In+l 



+ 



(t. 



n+1 



(3.2) 
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with Tqi := ti — to- Note that here the in (hence the Tnm) now play the role of parameters, since 
they are just the background solution. The actual perturbed time variables have disappeared 
from ()3.5p . mirroring the invariance of the background solution g„ = under translations of the 

in,' 

In the second line of ()3.5p we introduced a general parametrization of an action, quadratic in 
the variables xq,xi and symmetric under exchanging xq,xi. This allows also to consider different 
discretizations in one go. We will define a renormalization procedure (corresponding to the so- 
called decimation procedure, where one just integrates out some variables) on the space of such 
actions and look for fixed points of the renormalization flow acting on the parameter functions 
afc. Hence we can expect that some choices for the initial discretization are not relevant, as 
these different choices may flow to the same fixed point. 

A difference to the usual discussions of renormalization group transformations is, that the 
parameters afc(Toi) in the action (j3.5p are not just coupling constants but functions of the time 
distance Tqi. The appearance of 'coupling functions' a{Toi) can be avoided by expanding these 
functions in a suitable basis, for instance a power series in Tqi. In general this will introduce 
infinitely many coupling parameters (which one could also treat perturbatively). However the 
fixed point conditions (which are related to the condition of discretization independence and 
reparametrization invariance) will fix almost all of these coupling constants. 

To define the iteration procedure we assume that the one-step propagator has the form 



K{xo, xi,Toi) = r]{Toi) exp 



-- {ai{Toi){xl + xf) + q;2(Toi)xoXi) 



(3.6) 



that is that the measure factor t?(Toi) does only depend on the background time difference and 
not for instance on the xo,xi. This assumption is justified, as the form of the propagator will 
be stable under iteration. 

Since the perturbative variable ri does not appear in the action (|3.5p the propagator for two 
time steps 



K^^\xo,X2,Toi+Tu) = J dxidn ??^°^(Toi)r?(°)(ro2 

1 



exp 



5(0) (xi ,Tl,X2,T2) 



(3.7) 



leads to a divergent result. 

This infinite factor represents the volume of the discrete diffeomorphism group (in this case 
the actual placement of the interior discretisation point), which is the whole real line. We will 
drop this infinite factor here, so we will not integrate over ri, which will lead to a finite result^ 

Having dropped the ri-integration in the path integral, the remainder can be easily com- 



*This is reminiscent of the Ponzano-Regge model, which also needs to be gauge fixed in order to produce finite 
results. One could argue that ri shouldn't range over the whole real line, since ti should always remain between 
to and t2. fndeed, restricting the integration range appropriately would lead to a finite result, just as a similar 
restriction within the Ponzano Regge model - which ensures that one only sums over diffeomorphisms which 
preserve the orientation everywhere - renders the model finite. However, the finite result does not only not agree 
with the Ponzano-Regge amplitude, it also leads to a triangulation-dependent model, in which the diffeomorphism 
symmetry is still broken [57] • We therefore do not adapt that strategy here. 
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puted. Choosing T := Tqi = T12 for simplicitjH, we obtain for the propagator 



iC("+i)(xo, , X2, 2T) = r^^^\Tf j dxi exp (ai"^(r)(x2 + 2x1 + 4) + 4"^(r)(xoXi + X1X2] 



.^('^+i)(2T) exp 



( (2T) {xt + xi) + (2T) xqx^ 



(3.8) 



This defines recursion relations for the coefficients ai, 02 and rj given by 



(n), 



8aS")(r) 



(n+l) 



(2T) 



(2r) 



4"^(T)^ 
4a[")(T) 

2aS")(r) 



7?(")(r)^ 



(3.9) 
(3.10) 
(3.11) 



To find the perfect propagator, one can iterate the equations ()3.9p - ()3.1ip . with initial values for 
the af*^ (T) taken from (j3.5p . Alternatively we can directly look for the fixed points. Considering 
the first two equations involving only the a coefficients, a family of fixed points is given 



S* {xo,xi, T):=al(T){xl + xf) + a2(T)xoXi 
uj cosh(£l'T)(xo + x\) — 2xqXi 



2g 



sinh(a'T) 



(3.12) 



Note that this action is Hamilton's principal function for the harmonic oscillator with frequency 
oj and a coupling 1/g in front of the action. These constants are determined by the initial values 
for the action, which for our choice (j3.5p leads to oj = co, g = 1. Indeed [28], this initial action 
converges to ()3.12p under the iterations defined by ()3.9p . That there is at least a two-parameter 
family of fixed points parametrized by a), g can be easily deduced from the iteration equations 
p. op . These are invariant under a rescaling Ok — > g~^ak- Furthermore if ak{-) is a fixed point, 
so is afc(tD X •). In the appendix|B]we will show that the action ()3.12p provides the most general 
solution to the fixed point conditions. 

For the fixed point equation of the measure factor rj{T) we can use the fixed point solution 
aJ(T). A solution is given by 



r?*(r) 



2iThg sinh(tDT) 



exp(|r) 



(3.13) 



where ^ is a free parameter. Again we will show in the appendix iBl that this is the most general 
fixed point solution. The free parameter ^ describes the possibility to add a constant potential 
term ~ T to the action p.Sp . 

Note that the iteration equation for the measure factor r] is not linear in r/ - rather we have 
r] quadratically appearing on the RHS of (|3.1ip . Hence if we scale the initial value rj^^^ with a 



®It will turn out that this regular subdivision is sufficient to regain the full symmetry under vertex translations 
and arbitrary subdivisions, i.e. discretization independence. 

®The fixed point property is easy to check using sinh(2j/) — 2 sinh(i/) cosh(y) and cosh(2jy) = cosh^(j/)+sinh^(j/). 
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factor a we obtain a scaling for the n-th iteration. Therefore many initial values will 

either diverge or converge to zero. Starting with e.g. the measure factor of the path integral for 
the free particle r]^^\T) = \J T~2 leads to a non-vanishing convergent result however [28]. 
This also satisfies the normalization condition mentioned in section [2j 



To summarize for the initial action (j3.5p we obtain the perfect quantum propagator 

io cosh(wr) {xl + xf ) - 2xoXi ] 



K{xo,xi,T) = J ^ , \, exp 
" 27rnsmh(u;i ) 



h 2sinh(a;T) 



Note that this propagator (|3.14|) - although calculated for the linearized theory - coincides 
with the non-perturbative quantum mechanical propagator for the harmonic oscillator 



LP cosh(a;(ti - to))(go + gl) - 2gogi 
h 2 sinh(a;(ti — to)) 



(3.15) 



This also gives the propagator for the parametrized harmonic oscillator - which we obtain from 
our calculation for the linearized theory by replacing the background parameters tn by the full 
dynamical variables t„ and the perturbations Xn by qn- 
Hence the propagator property 

K{qo,to,q2,t2) = j dqi K{qo,to,qi,ti) K{qi,ti,q2,t2) (3.16) 

is not only satisfied for equal step times {t2 — ti) = (ti — to) (which holds due to the fixed point 
condition) but also for arbitrary subdivisions into unequal time steps0 Furthermore we achieved 
full 'discretization independence': if we use for the discrete propagator in the (discretized) path 
integral the perfect propagator the result will not depend on the number of subdivisions we use. 

As explained in section [21 related to this discretization independence is the invariance of the 
path integral (j3.16p under 'vertex translations', i.e. changing the value for ti in (j3.16|M Again, 
we do not integrate over ti (which in the parametrized theory is a dynamical variable) in (j3.16p . 
as this would just lead to an infinite factor. 

The gauge symmetry in the path integral leads to a constraint in the canonical theory. Indeed, 
the propagator ()3.15p satisfies the quantum mechanical constraint equation for the parametrized 
harmonic oscillator (which is just the usual Wick rotated Schrodinger equation if one reinterprets 
t as a non-dynamical time parameter) 



d d 



2 ^2 



CK{qo,to,qi,ti)= [h— + —-^ + —q'ojK{qo,to,qi,ti)=0 . (3.17) 

Note that this constraint equation is only satisfied for the fixed point ()3.15p and not for any 
of the propagators K^""^ for finite n, which would correspond to some (non-perfect) discretized 
version of the path integral. That is the fixed point is characterized by a symmetry, which is 
reflected in the constraint equation ()3.17p . 



'^This is not so surprising as we basically computed the path integral iteratively. The fixed points correspond 
to taking the (equal time step) discretization to the continuum limit. 

**The gauge symmetry also involves a corresponding change in qi, this can however be absorbed in a variable 
transformation for qi in the action. One can check that the Jacobian of this transformation changes the measure 
factor accordingly, so that one obtains the measure appropriate for the transformed ti. 
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3.1 Linearized dynamics around a solution with qn ^ 

In this section we will linearize the discretised, parametrised harmonic oscillator around a clas- 
sical solution ()3.3p with ijn 7^ 0. Since in this case the background solution does not exhibit 
gauge symmetry, neither does the linearized theory [H]. Hence, the appear as variables and 
have to be dealt with. 

Expanding the action (|3.'2|) into second order via (|3.3|) . one arrives at an action for the 
linearized variables Since reparametrization symmetry is broken, there is no obvious 

choice for a gauge variable (a role which was fulfilled by the r„ in the previous case, where the 
symmetries were present), so in the path integral the integrals over both Xn and t„ have to be 
performed. Moreover, the resulting propagator is not a projector on the physical Hilbert space 
of the continuum theory. 

Performing one refinement step therefore amounts to evaluating 

i^("+i)(a;o,ro,X2,r2) = / dxidrii^(")(xo, tq, xi, ri)i^(")(xi, n, X2, r2) (3.18) 



with K^^\xo,To,xi,Ti) = rj^^^ ex.p{—S^^\xo,To,xi,Ti)/h). In this case r/(") and the action 5^"^ 
depend furthermore on the background variables qi,ii. It should be noted that, for one refine- 
ment step, the background solution qi, ti has to be chosen such that it is a solution to the 
equations of motion of the discretized action 5^"^ with boundary values qo, to, q2, ^2- 

The integral (j3.18p is actually finite in most case^, despite the ri integration actually being 
carried out. This is a result of the broken symmetries of the linearized theory: the t„ are actually 
propagating degrees of freedom, corresponding to a small but non-vanishing eigenvalue in the 
Hessian of the action. 

Due to Ti being present, there are many more coefficients in the action S^^\ leading to more 
complicated recursion relations coming from ()3.18p . These have been derived in [28j, and we 
report on the findings in what follows: 

The recursion relations for the coefficients coming from the action 5*^"^ converge to a fixed 
point which is given by 

S = + ^) - + ,^^(co.h(.T)(4 + 4) - 2.,,.,) (3.19) 

2 smh(cji j 2smh(ci;jj ' 

w2(go - cosh(cjr)gi) a;2(gi - cosh(tjr)go) . ^ 

+ . ,2/ ^o(ro - Ti) + ■ .2/ ^1(^0 - Tl) 

smh [ujl ) smh [oji ) 

io^{qo — cosh{ujT)qi){qi — cosh(a;T)^o) , 2 2 



2sinh*^(a;r) 



T^ + rt) (3.20) 



with T = ti—tQ. Not surprisingly, the action (|3.20p is actually the perfect action of the harmonic 
oscillator expanded up to second order in the linearized variables ()3.3p . More interesting is the 
behavior of the measure factor ri^^\ One can separate the measure factor into two contributions 



^Note however that even in the Wick-rotated framework the path integral for the hnearized theory H3.18P is 
not necessarily always convergent. The reason is that, depending on the (background) boundary values go, ^o, 
Qjv, tpf, the Hessian matrix of second derivatives, which arises as inverse propagator for the x„, t„ due to the 
expansion H3.3p . might not be positive definite, rendering the Gaussian integral divergent. This is also an effect 
of the broken gauge symmetries. There is, however, a good control over which boundary values lead to positive 
definite Hessian [28], and we assume that such a choice has been made in what follows. 
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where A^"^^ arises as prefactor from the integration over the xi, and B^"^^ from the integration 
over Ti in ()3.18p . It is not hard to show that in the limit of infinite refinement 



^^"^ ~^ y 27Thsmh{ioT) ^^-^^^ 

which coincides with the measure factor for the linearized case (j3.14p . However, despite S^") 
being finite for every n, in the limit it actually diverges, leading to a diverging propagator. 

This is also not surprising if one keeps in mind that in the continuum theory r is actually the 
gauge variable, and therefore should not be integrated over in the path integral. The breaking 
of symmetry in the discrete theory has led to a finite integral, however in the perfect limit the 
integration over ri coincides with an integral over the volume of the gauge group. This volume 
is infinite, since it corresponds to the placement of the intermediate point t„,, which can be put 
everywhere on the real line. 

Hence, in the perfect limit, the gauge symmetry can also be recognized by the volume of 
the gauge group still being present in the path integral. This could have been remedied from 
the first place by not integrating over ri, but rather gauge fixing it to a definite value, e.g. 
Ti = {tq + T2)/2, which also coincides with the classical solution to the equations of motion of 
the linearized system. In the case of gravity such a split is not so easy to obtain. In this case one 
has to consider the eigenmodes of the Hessian of the action and the corresponding eigenvalues. 
The eigenvalues of the pseudo gauge modes will converge to zero in the perfect limit and will 
lead to a divergence of the measure. 

Dropping the infinite volume of the gauge group, one arrives at the perfect propagator 



K{xo,ro,x„n) = y^^^^J^|^exp(-i5(xo,ro,xi,n)j (3.23) 

where S is given by ()3.20p . Is is straightforward to check that the propagator ()3.23p satisfies the 
constraint 

C K{xo,TQ,xi,Ti) = [h- h —^r-^ + —xl] K{xQ,TQ,xi,Ti) = (3.24) 



dn 2 dx{ 2 

i.e. K \s a. projector on the physical Hilbert space. Furthermore, writing explicitly the depen- 
dence of K on the background solution qi and tj, K also satisfies 

K{xq, To, X2, r2; to, q2M) (3.25) 
= j dxi K{xo,To,xi,Ti;qo,to,qi,ti)K{xi,Ti,X2,T2;qi,ti,q2,t2) 
for every qi, ti that satisfy 

_ sinh(t2 -h) , _ sinh(ti - to) „„x 
smh(t2 - to) smh(t2 - ^o) 

i.e. that satisfy the equations of motion of the (Wick-rotated) continuum harmonic oscillator. 
In this sense the perfect propagator is independent of the actual discretization, i.e. for N = 2 
of which actual background solution one is perturbing around. 
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4 The anharmonic oscillator 



In section [3] we dealt with the harmonic oscillator linearized around the background solution 
with g„ = 0. Here the linearization effectively lead to omitting the fluctuations r„ in the time 
variables. Indeed, treating these explicitly in the path integral is quite complicated, as the time 
variables appear non-polynomially in the action. 

To discuss a truly interacting theory we will now consider the quartic anharmonic oscillator, 
so that we need to expand the action to at least fourth order. We have to face the problem 
that the r„ appear in the third and fourth order of the expansion. Nevertheless we will see that 
also in this case we can avoid the integration over the fluctuations r„ by perturbing around the 
perfect propagator for the harmonic oscillator. 

To start with, a discretization for the quartic anharmonic oscillator is given by the following 
action for one time step 

Soi = l^-^^ + Y^ll + ^o)(*i - *o) + ^(^0^ + - to) . (4.1) 

In a path integral approach we can treat such actions only perturbatively. For instance we can 
expand as before around the background solution = and i„ arbitrary. To second order in 
the perturbation variables Xn,Tn we recover the harmonic oscillator, treated in section [3l In 
particular there the perturbations did not appear and we could just ignore the integration 
over these variables in the path integral. 

To capture however the anharmonic term, we need to expand to higher order in the variables, 
which introduces a dependence on the Tn- Even worse, the lowest order terms in which the r„ 
appear, are linear in the t„, which cannot be directly dealt with in the path integral. 

Classically one can show |14] that a perturbative treatment of actions with broken sym- 
metries, which are however expanded around a background solution with exact symmetries, is 
in general inconsistent in higher-than-linear order. The problem is that the background gauge 
parameters (here the tn) are free to lowest order in perturbation theory but are fixed by (non- 
linear) consistency conditions to higher order. That is, if the t„ are not already chosen such 
that these consistency equations are satisfied, the ansatz t„ = tn + eT„ is not justified as the 
solutions for the Tn will involve terms Tn ~ e~^. 

To avoid this issue one can improve the action order by order. That is one would start with 
the quadratic order of the action and find the perfect action to this order. To this one adds the 
(naively discretized) third order term. For this third order action the problem above does not 
appear (to lowest non-linear order), i.e. the background gauge parameters i„ remain free. This 
third order action has to be improved again, after which one can add the fourth order term and 
so on. 

This solves the consistency problem on the classical level. It also solves the problem in the 
path integral: the action (and measure) being perfect up to a certain order will be also gauge 
invariant to this order. Hence one can either apply gauge fixing, or rewrite the action into gauge 
invariant variables, such that the path integral needs only be performed over these variables. 

As explained in section [3] in the case of the harmonic oscillator the perfect action for the 
quadratic approximation gives us already the full non-perturbative perfect action, if we just 
replace the background time parameters i„ with the full variables tn- We therefore do not need 
to perform the perfection procedure to third order. To fourth order we have a contribution from 
the anharmonic potential term. So we will start from the action 

a _^ f cosh(a;(ti - to))(gg + gf) - 2qoqi \ A . 4 , 4w , ,^ f.^-. 
- 2 [ sinh(a;(ti - to)) j + 2^4! + ^^^^^^ " ' ^^'^^ 
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We will again consider the case with two time steps, so that we always need to integrate only 
over one variable pair qi,ti. That is we consider 



S = Soi{qo,to,qi,ti) + Si2{qi,ti,q2,t2) ■ (4.3) 

Applying the expansion = + 3;„ , t„ = t„ + t„ we would actually encounter third and fourth 
order terms (from the expansion of the perfect part of the action) in which the ri appear. We 
know however that the perfect action is exactly gauge invariant, hence there exist a variable 
transformation ((?i,ti) (Qi,Ti) such that this perfect part only depends on Qi. Indeed with 
Ti = ti and 



/ sinh(a;(t2 - tp)) _ sinh(u;(t2 ~ ^i)) ^^ _ sinh(^(ti - ^o)) ^^ i 

y sinh(cj(ti — to)) sinh(a;(t2 — ti)) V sinh(a;(t2 — to)) sinh(cj(t2 — to)) 



we find 



c , c , w /cosh(a;(t2 - to))(go + ^i) - 2gog2 

JQl + '-'12 — TT*^! + 



2^^ 2 V sinh(a;(t2 - to)) 

+ ^Uqo+QfiQu<lo,q2))ih-to) + {qf{Qi,qo,q2)+ql){t2-h)) . (4.5) 

In the last line qi has to be expressed as a linear combination of Qi^qo and q2 (with t^ 
dependent coefficients) by inverting ()4.4p . Note that the second term in the action (|4.5p is 
just Hamilton's principal function for the harmonic oscillator, coinciding with the perfect action 
for the time step t2 — ^o- We now expand the action (j4.5p in Qi = + Xi, g„ = + x„ for 
n = 0,2 and in t„ = t„ + r„. Indeed we see that ri does not appear in an expansion up to 
fourth order in the variables. We can therefore drop the ri-integration in the path integral just 
as in the case of the harmonic oscillator. A similar argument can be applied if we consider the 
path integral involving any number of steps. Again one can find a transformation from the 
qi,. . . ,g„_i,ti, . . . ,tn-i to variables Qi, . . . ,QN-i,ti, . . . jt^v-i such that the harmonic part of 
the action depends only on the Qk- We can therefore ignore the fluctuation variables also in 
the iteration process (which computes the A^ = 2^^^ path integral iteratively) below 

This allows us to compute iteratively the perfect propagator for the anharmonic oscillator to 
first order in A. As for the harmonic case we have to choose a parametrization for the propagator. 
To this end consider the first iteration step, i.e. integrating out Qi from (|4.5p . where again we 
assume T = t2 — ti = ti — to: 

K(^\xo,X2,2T) = j dQi{7]{T)f exp(^-^Ql-^Sharm{xo,X2,2T)^ X 

;-^(4?(T) + /i?(T))+0(A^)) 



:r?(2T)exp ( -^Sharm{xo, X2,2Tyj (^1 - ^lg\2T) + 0{X^) ] (4.6) 

The action Sharm is the perfect action for the harmonic oscillator appearing in (|4.5p . It does 
not depend on Qi, hence this exponential factor can be pulled out of the integral. 



Alternatively [2H1 to introducing the variable Qi one can expand first the action 
for n = 0, 1, 2 to fourth order and then define a coordinate transformation xi Xi such that the action to fourth 
order only depends on Xi and not on ri. This leads to the same results as the approach presented here. 
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We have chosen a measure factor 



vm= I " ( ""^^^"^ (4 7) 

which corresponds to the perfect measure for the harmonic oscillator for the Qi variable. The 
second factor is (the square root of) the Jacobian of the transformation from qi to Qi. To zeroth 
order in A this measure factor remains indeed invariant under iteration. (First order corrections 
in A to the measure can be absorbed into the interaction term /o2-) 

The interaction terms , are polynomials of up to fourth order in Xi. The coefficients 
of the Xi in these polynomials depend on xo,X2 in such a way that all terms are fourth order 
if we add up the powers in xq, Xi, and X2. Performing the integratior0 in ()4.6p one will find 
that /q2'' '^ill contain all even powers in xq, X2 up to fourth order including a constant term. We 
therefore adopt the following parametrization for the interaction term 

Ioi=ao{T) (xq + xf) + ai(r) (xgXi + xoxf) + a2{T)xlx1 + 

/3o(T) {xl + xj) + /3i(r) xoxi + 7(r) . (4.8) 

From the iteration 

K("+i)(xo,X2,2r) = j dXiK^'^\x^,xi{Xi),T) K^''\xi{Xi),X2,T) + 0{\^) (4.9) 

with 

K(")(xo,xi,r) = ry(r)exp (-^^^^(xo, xi, T)^ {l - ^I^^ {T) + 0{\^)^ (4.10) 

we obtain recursion relations for the coefficients a, /3, and 7. This is a straightforward exercise 
involving only the transformation from xi to Xi given by ()4.4p . just replacing qi, Qi by xi, Xi 
and performing the integration over Xi. 

By expanding in powers of h these recursion relations can be divided into a classical part and 
quantum corrections. The classical part of the recursion relations coincides with the recursion 
relations one would obtain by extremizing 

Sharm{xo,X2,2T) + XI^'^^'^^Xq, X2,2T) = Sharm{xo, Xi,T) + Sharm{xi, X2,T) + 

A/(")(xo, XI, T) + A/(")(xi, X2, T) + 0(A2) (4.11) 

with respect to xi. A solution to the corresponding fixed point equations is given by Hamilton's 
principal function (for the corresponding continuum action), which can be readily obtained to 
first order in A, see appendix IK\ 

We will see that the classical part coincides with the homogeneous part of the recursion 
relations. As will be explained in the appendix ^ the inhomogeneous recursion relations can 
be easily brought into a standard form once the fixed point solutions to the homogeneous part 
is known. Hence it helps very much to have the classical system solved, in order to obtain the 
full quantum mechanical solution. 
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The recursion relations for the aj-coefficients are a closed system as these correspond to the 
classical problem with just an order interaction term: 

4"+^)(2r) = ^ . 4 . f (1 + 8cosh^(ru;)) 4")(T) + (4cosh3(r..) + cosh(T.;)) ai")(r) + 
8 cosh [Tlo) V 

2cosh2(rw)4"^^ 

S"^'^(2T) =---l^(^4")(r) + (cosh3(r.;) +cosh(T^))aS")(r) +cosh2(ru;)4'^)) 

(2T) = ^^^2\Tco) (34"^m + 3^°^h(T^)ai")(T) + 2cosh^(Ta;)4")(T)) . (4.12) 



a 



(n+l) 

On 



al(T) = — r-TTT^i ~ ^'^Tuj cosh{Tu) + 9sinh(ra;) + sinh(3Tw; 



One fixed point of these equations is provided by the perfect action for the anharmonic 
oscillator (to linear order in A), determined in the appendix lAl 

an(T) = fl2Ta; - 8sinh(2ra;) + sinh(4Tw^ 

°^ ^ 768cjsinh^(rw) V ^ ^ ^ ' 

192a; sinh*(ra; 

a*(T) = ^ f2ra;(2 + cosh(2ra;)) -3sinh(2Ta;)') (4.13) 

^ 64a;sinh^(TL<j) V ^ ^ ^ V ' 

The fixed point equations are invariant under a rescaling — > Xoi, and hence we have A 
as a free parameter for the solutions. This is easy to understand as the freedom to rescale 
the interaction term by an arbitrary constant, redefining the coupling constant A. As for the 
ambiguities g, Co appearing in the fixed points for the harmonic oscillator, the final coupling 
constant is determined by the initial conditions for the iteration procedure. 

The solution ()4.13p is actually not the most general one. We will explain in the appendix [B] 
that the most general solution has two further free parameters, which determine the couplings to 
terms and iP'x^ in the corresponding continuum Lagrangian. In the following we will however 
set these couplings to zero, that is consider the standard quartic anharmonic oscillator with a 
perturbation term x^. 

Next we consider the recursion relations for the /3j, which are given by 



,i"«>(2T) = iT) (l + ^^^) + (T) 

^(n+i)^2T) = /3(") (T) + /3i") (T) ^ 



cosh^(ra;) ' cosh(rw) 
^ . 3tanh(ra;) („) 3tanh(ra;) („) x 

oj\cos\i^{Tuj) ° ^ ' 2cosh(Ta;) ^ ^ V ^ ^ 

The equations ()4.14p have a homogenous part which is independent of ^, and an inhomogeneity 
proportional to h depending on the Oj, i.e. the fourth order potential. To any solution of the full 
recursion relations one can add an arbitrary multiple of the solutions to the homogeneous part 
of the recursion relations. Hence to find all solutions one has to determine also the solutions to 
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the homogeneous part. As this part is of the order of it is again equivalent to the classical 
iteration relations one would obtain if one considers an interaction term quadratic in x added 
to the perfect action for the harmonic oscillator. A solution for the homogeneous fixed point 
equations is hence again given by the perfect action for the harmonic oscillator: Take (j3.19p and 
replace there oj, g uj + XD, 1 + A/i and expand to first order in A, resulting in 

/3^(T) = ^coth(Tw) ''^'^ 



sinh^(Ta;) 

PHt) = V 0, \ - M — T77^ ■ 4.15) 

' smh^{Tuj) ^sinh(rt^) ^ ' 

(Here we reparametrized /i = — w/i) and v = ^P.) Again we find a two-parameter ambiguity 
corresponding to the ambiguities g and w that we found for the harmonic oscillator. 

For the inhomogeneous equations we can assume that the ai are given by their fixed point 
values (j4.13|) . There is a general strategy with which one can attempt to find the fixed points 
for the inhomogeneous equations, which requires the knowledge of the homogeneous equations. 
We will demonstrate this in the appendix O as an example. The general solution is given by 

/3o*(T) = /icoth(ra;) - -^^^^ 



sinh (To;) 



1^ .2/ 



+ ^ 97 r 2 + cosh^fTcj) - STo; cothfTo;) 

32a;2sinh2(ra;) V ^ ' ^ ' 

2Tujcosh{Tuj) 2 
= u — . ^ ^ fi 



smh^iTu) '^sinh(rw) 

H f ^Tcj + 2Tuj cosh(2rw) - 3 sinh(2Tw)) . (4.16) 

32a;2 sinh (Tw) V / 

Finally the recursion relations for 7 are given by 

^(n+i)(2T) = 27(")(T) + ^tanh(ra;) ^(„)^^^ ^ 3^W(r^^H^^^ _ 



Again we have a homogeneous part corresponding to a classical iteration procedure. This corre- 
sponds to having a constant potential added to the perfect action for the harmonic oscillator. It 
is easy to see that under the classical recursion relations such a constant term is just multiplied 
by 2 corresponding to changing the time interval from T to 2T. Hence a (family of) fixed points 
for the homoegenous part of the relations is just given by ^*{T)hom = CT. 

The fixed points for j*{T) for the full recursion relations can be readily found with methods 
shown in appendix O These are given by 

r(r)=^r-^(;.-.r.coth(ra;))-^(3coth(rc.)-Tc.(2 + -^^ . (4.18) 
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To summarize, the fixed point propagator is 



27rnsmh[I uj) \ n J \ a 

A 



768a; sinh"^ {Tuj 
A 

192w sinh^(ra;) 



l2Tuj - 8sinh(2ra;) + sinh(4ra;)j (xq + x\) + 
- 12ra;cosh(ra;) + 9sinh(Ta;) + sinh(3Ta;)^ (xqXi + xoxf)+ 



^ ^ ^ (2Tu;{2 + cosh(2Tw)) - 3 sinh(2ra;)) xlxl+ 



64a;sinh (Tw) 

j^j^^ (2 + cosh^ (Tuj) - 3Tuj coth(Ta;) 

u coth(Tuj) ^- + \h ^- ^ ) (xn + 2;?)+ 

^ ^ ' sinh2(rw) 32a;2sinh2(ra;) ° ^' 

2Tuj cosh(Tuj) 2 ~ (4r6j + 2Ta;cosh(2rw) -3sinh(2rL^) 

y /i -, r + A/l-^ 5 — XnX^ + 

sinh2(rL^) ^sinh(rw) 32w2sinh3(rtj) i 1^ 



iT --(11- vTuj coth(Tw)) - (3 coth(Tw) - Tuj(2 ^ 



UJ 



64a;H ^ ' ^ sinh2(ra;)' 



+ 0(A 



2\ 



(4.19) 



For the initial data used in (|4.6p the propagator takes on a simplified form, as in this case 
A = 1 and n = v = = 0. Nevertheless the result is quite complicated: even to guess the correct 
reparametrization invariant measure (i.e. all terms proportional to h and h"^ in square brackets) 
in case the perfect classical action is given, without actually solving the dynamics, seems to be 
quite impossible in general situations. 

From the infinitely many parameters ai{T), /3i{T), and j{T) in our initial parametrization 
we are left with the couplings A, u, and ^ (and a further two couplings corresponding to 
adding terms and x^x^ to the continuum Lagrangian, see appendix [B|) . These parameters 
characterize the continuum Lagrangian (and Hamiltonian, see below), i.e. all discretization am- 
biguities are resolved by requiring reparametrization invariance for the discretized path integral. 



Again, although for the recursion equations (|4.9p we have set ti = the fixed points of 

the recursion equations satisfy a stronger condition, namely 

K{xQ,X2,Ti+T2) = j dxiK{xo,xi,Ti)K{xi,X2,T2) + ©(A^) (4.20) 

for arbitrary positive Ti, T2. In other words, the perfect propagator K{xo, xi,T) leads to a 
path integral which for > 1 has become independent - to first order in A - on the actual 
placement of the intermediate discretization points, which correspond to the discretized choice 
of t{s) in ()3.ip . It therefore mimics exactly the gauge symmetry of the continuum theory. Also 
()4.20p shows that the model defined by the propagator ()4.19p is discretization independent (up 
to terms of order A^), i.e. the discrete path integral defined by the amplitude (j4.19p does not 
depend on the number of discretization points. 
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The invariance under a gauge symmetry leads to a constraint, which is satisfied by the 
propagator (|4.19p 

,^ h-^ + {l + 2{v - M))y £2 + i^^^ + + ^)) + ir^i + ^ (^-^^^ 

such that 

C^'^'^'^'^^i^lxcXi.T) = 0(A2) . (4.22) 

This constraint equation characterizes the fixed points, it would not be satisfied by the K^"'^ for 
finite n. 



5 Summary and Discussion 

In order to obtain a well defined path integral one often utilizes discretizations. This may 
however lead to a breaking of (gauge) symmetries, which in the case of general relativity are 
central for the dynamics of the system. We addressed this problem here in the context of 
reparametrization invariant systems. Despite being extremely simple compared to gravity, these 
systems share the property that the gauge symmetry determines the dynamics of the system 

m- 

Indeed we could show in section [2] that requiring the implementation of reparametrization 
invariance for the discretized path integral - where it assumes the form of vertex translations 
(|2.4p - uniquely fixes the discrete propagator. Namely the propagator has to coincide with the 
quantum mechanical continuum propagator of the system under consideration. Furthermore 
such a discretized path integral is automatically discretization independent, i.e. the propagator 
does not depend on the number of subdivisions . This is just the convolution property (j2.7p of the 
continuum quantum mechanical propagator. The number of subdivisions can be taken to zero, 
in which case one obtains the discrete propagator itself, even for large time steps. This (fixed 
point) discrete propagator is furthermore characterized by satisfying the constraint equations. 

We conjecture that similar properties will also hold for discrete gravity: a full implementation 
of diffeomorphism invariance in the form of a symmetry under vertex translations [22J should 
lead to discretization independence. Similar to the example of reparametrization invariance, the 
discretization can then be taken to be a very coarse one, hence one would expect that many 
discretization ambiguities will be fixed (as free parameters in the discrete propagator have an 
effect on macroscopic scales). It would be interesting to have a proof for more general cases, 
in analogy to the arguments presented here. One important difference to the one-dimensional 
case is that a perfect action or perfect discretization for higher dimensional theories will include 
non-local couplings |2n1 [2T| [22] . 

The implementation of vertex translation symmetry into the discretized path integral will also 
ensure that this path integral satisfies constraints [13], which can be taken as a characterization 
of the fixed point. Again, an important difference to the one-dimensional case is, that such 
constraints will involve non-local couplings, and are therefore not explicitly known. These 
constraints would however be free of anomalies. 

Coming back to reparametrization invariant systems, we can take the requirement of dis- 
cretization independence in the form of the convolution property ()2.7p as a starting point to 
define an iterative procedure. This iterative procedure improves the discrete propagator, so 
that in the limit it satisfies the convolution property (|2.7p and defines a discretization in which 
reparametrization invariance is realized in the form of vertex translations. Some points we wish 
to emphasize, are the following: 
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Quite naturally, one is led to consider not one specific discretization, but an entire class 
of discrete models. The iterative procedure defines a renormalization group flow within 
this class. This allows a discussion of all possible fixed points. Indeed here, we found not 
only the solutions to the (an-) harmonic oscillator, but to the most general Lagrangian 
involving even powers in q and q up to order 4. Hence to study the relevance of discretiza- 
tion ambiguities in gravity models |TT], it might be essential to study the behavior of these 
models under coarse graining. Note also that the relevant parameters - which fix the corre- 
sponding continuum Lagrangian - are all determined by the lowest orders in an expansion 
of the discrete propagator in the time variable, that is by the behavior of the propagator 
for short times. The higher order coefficients in this expansion can indeed be understood 
as proper discretization ambiguities, which are fixed by requiring reparametrization invari- 
ance of the discrete propagator. This might be an interesting point for spin foam models, 
as there one often rather concentrates on the limit of having large building blocks (with 
an amplitude that would correspond to what we termed naive discretization, as it has not 
been subject to any coarse graining procedure). 

Even for the simple case of the anharmonic oscillator, the reparametrization invariant 
discrete path integral (|4.19|) is quite complicated, and it would probably be impossible to 
guess it without actually solving the dynamics. Similarly, to find an anomaly-free measure 
for spin foam models [HI [171 [H]) it seems to be unavoidable to address the dynamics of 
the system, in particular the behavior under coarse graining. Even to just ensure anomaly- 
freeness in the continuum limit (which is the only thing one might realistically expect), a 
study of the coarse graining properties might be valuable. As we have seen, the choice of 
initial data for the iterative procedure of the measure factor had to be done carefully in 
order to obtain a convergent result. The behavior of the amplitudes under (dynamically) 
trivial subdivisions |161 [18] are nevertheless interesting as first steps in this direction. These 
might actually be important in order to show that vertex translation symmetry implies 
discretization independence for gravity, in analogy to the one-dimensional case. 

Path integrals with (properly implemented) gauge symmetries act as projectors onto the 
space of physical states [13]. Indeed spin foams are often mentioned as a tool to obtain 
physical states, or to define the physical inner product, for loop quantum gravity. As we 
have seen in the toy example of reparametrization invariant systems, to obtain a propagator 
satisfying the quantum constraints, it was necessary to take the fixed point propagator, 
that is the perfect or continuum limit of the (naively) discretized path integral. On the 
other hand the constraints can be expressed as conditions on the fixed point coefficients 
a|(T), . . . ,7*(r). For future work it would be interesting to explore the relation between 
these conditions encoding reparametrization invariance, and the fixed point conditions 
encoding discretization independence. 

Finally, the methods applied in this work - basically a version of Wilsonian Renormaliza- 
tion Group flow - might be actually useful in order to find solutions to quantum mechanical 
path integrals. Also it would be interesting to apply these techniques to reparametrization 
invariant systems arising in mini-superspace reductions of gravity, for instance for loop 
quantum cosmology |3U] . The path integral for the anharmonic oscillator is usually used 
to derive the corrections to the energy levels of the harmonic oscillator ^J, for which 
one just needs to obtain the T ^ oo behavior of the propagator. Here we derived the 
full propagator for arbitrary T, as this is needed to define a perfect discretization of the 
path integral, and with this an anomaly-free (with respect to vertex translation symmetry) 
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path integral measure. We also want to point out the work [32], where the authors also 
discussed a perfect path integral for the anharmonic oscillator. The difference to the work 
presented here is that [32j uses a different coarse graining procedure, namely averaging 
instead of decimation as applied here. We employed decimation as this seems to be the 
only method to obtain reparametrization invariance. 
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A Hamilton's principal function for the anharmonic oscillator 

Here we give Hamilton's principle function for the anharmonic oscillator to first order in A, as 
this will provide a fixed point for the recursion relations for Oi (j4.12p . 

In [3J it was shown that for a ID system the perfect action coincides with Hamilton's principal 
function for the given boundary values, i.e. 

Sperf{qo,to,qi,h) = J ds (^^^ + ^qH' + ^qH'^ (A.l) 

where q{s) and t{s) are solutions to the continuum equations of motion with boundary values 
QO: to, Qi, and ti. To find Hamilton's principal function to first order in A we would need to find 
the solutions at most to first order in A, and then perform the integral in (jA.ip . i.e. we expand 

q{s) = qis) + Xxis) + 0{X^) (A.2) 
t{s) = iis) + Xt{s) + OiX^) (A.3) 

where q, t are solutions to the harmonic oscillator. However to first order in A we do not need 
the explicit form of the solutions x{s), t(s): These would only appear in the harmonic oscillator 
part of the action (jA.ip . This contribution vanishes however due to the (harmonic oscillator) 
equations of motion for the background solution q{s), i{s). Hence we just need to evaluate the 
integral (jA.ip on the harmonic oscillator solution. One finds 

c- / . u} cosh{Tuj){q^ + qf) - 2qoqi 

^perf{qO,to,qi,ti — — + 

2 smh(j uj) 



A 



(^12Tuj - 8 sinh(2ra;) + smli{4Tuj)^ {q^ + qf) + 



768a; sinh*(rc< 

- ASTuj cosh{Tu}) + 36sinh(ru;) + 4sinh(3Tw) ) {qoqf + q^qi) + 



(24Tu}{2 + cosh(2Ta;)) - 36sinh(2Ta;)^ q'^qj 



(A.4) 



Similarly one can obtain Hamilton's principal function to first order in A for the harmonic 
oscillator with perturbation terms Xq^ and q'^q'^. As explained in the next section, these terms 
arise in the most general fixed point solution to the recursion relations (|4.12p . 
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B On the uniqueness of the fixed point solutions 



Here we will discuss the uniqueness of the solutions to the fixed point equations (j3.9p . (jS.lOp . 
()3.1ip . ()4.12p in the main text. In all cases we assume that the the solutions can be represented 
by a power series 'Yl'n=no ^nT"", which starts with some finite lowest power T"", that can also be 
negative. 

We start our considerations with the relations (|3.9p - (|3.1ip for the discretized action of the 
harmonic oscillator, which we rewrite into 

= al{2T) al{T) - a\{Tf + \al{Tf (B.l) 
Q = al{2T)a\{T) + \al{Tf . (B.2) 

Making the ansatz 



oo 
n=no 



(r)=X;a,,„r'^ (B.3) 



one will find the following equation arising from the coefficients to the lowest power T in 
the equations (jB.ip 

0= (2'^o-l)(ai,„„)2+ 1(^2, 



8 y^I^no) 

0= 2^^° ai^no a2,no + \{a2,n,f ■ (B.4) 



It is easy to see that this equation can be only solved for no = —1, in which case we obtain 

02,-1 = Ki , ai-i = (B.5) 

where ki is a free parameter. After having fixed the lowest order, one can convince oneself that 
if one iteratively solves the higher order equations for the coefficients of j'-^+fc+i^ then these 
equations are linear (inhomogeneous) equations for the coefficients aj^^, i.e. are of the form 



Here the matrix A is given by 



'^Aijaj^k = hi . (B.6) 



nfc-l _ 3 _i 
A={ _i ^ ) (B.7) 

2 ^ 2, 



and hi represents the inhomogeneous terms. The matrix has only vanishing determinant for 
k = 1 (and k = —1, which we already discussed). Indeed, the equations for aj,i add a further 
free parameter K2, as one will find the solutions 

"2,1 = K2 , = -K2 . (B.8) 

Since the matrix A has non- vanishing determinant for all other /c's one will have unique 
solutions for all the other coefficients depending on the two free parameters ki, K2- The solution 
obtained in this way agrees with the one (j3.12p presented in the main text, with an appropriate 
choice of the two free parameters g, uj there. 

One will find the same matrix A appearing in the recursion relations (j4.14p for the /3- 
coefficients in the anharmonic oscillator case. This is not surprising, as these arise by adding a 
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perturbation to the harmonic oscillator quadratic in the variables xq, xi. Since the determinant 
of the matrix A is only vanishing in two cases k = —1, 1 one will again find at most solutions 
with two free parameters, which is indeed the case for the solutions ()4.16p . 



We will now discuss the solutions to the fixed point equations (|3.11|) 

oo 



ri*{2T) 



1 7r?isinh(a)T) 



2 cosh(6DT) 



k=0 



(B.9) 



Here we used in the fixed point solution ()3.1f p and expanded the prefactor appearing in the 
first equation into a power series. As we have to choose for our initial measure r]^^^ a functional 
dependence rj^^^ ^ j'-i/2 ^-^g start with the assumption that ry* is of the form 



(B.IO) 



n=no 



with some finite, not necessarily positive number uq. Using this form in the fixed point conditions 
(|B.9p . it is easy to see that we must have uq = and that the first coefficient r]o is given by 



(B.ll) 



Using (|B.lip . one will find that the coefficient equation for the power of T^-i/^ leaves 
r]i as a free parameter. Indeed for the coefficients % with k > 0, the equations are of the form 



(2'' -2)rjk = fk{m,l<k) , 



(B.12) 



with /i = 0. Hence rji remains a free parameter, determining all other coefficients uniquely. 
The additional parameter corresponds to adding a constant potential term to the Lagrangian 
for the harmonic oscillator. The full solution is given in (|3.13p . 



Now we will turn to the recursion relations (|4.12p 
1 



«o(2r) 



al{2T) 
al{2T) 



. w,^ ^( (l + 8cosh^(r^))aS(r) + (4cosh3(r^) + cosh(ra;))at(T) + 
8 cosh (Tw) 



1 



2 cosh'' (To; ) 
1 

4cosh^(TLt; 



2 cosh^(rti;) a2 
al{T) + (cosh3(ru;) + cosh(ru;)) al{T) + cosh^{Tuj) a\ 
(3aS(r) + 3cosh(rw)a];(r) + 2cosh2(ra;)a^(r)) 



(B.13) 



for the terms in the discretized action of the anharmonic oscillator. This case is even easier 
to treat than the relations (jB.ip as we have a system linear in the variables. We assume again 
an ansatz 



a. 



(B.14) 



7i=no 
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with a lowest order T"". Since the cosh(Ttj) ' functions appearing in (|B.13p can be expanded 
in a Taylor series (i.e. there are no negative powers of T appearing), the equation for the lowest 



order coefficients ai^no will be 



(B.15) 



where 



A 



1 - 2" 

3 
4 



(B.16) 



The determinant of this matrix is only vanishing for k = —3, —1 and k = 1 and for these cases 
the matrix has rank 2. Hence we can expect three linearly independent solutions. For the higher 
order coefficients one has the same equation as in ()B.15p . just that inhomogeneous terms (arising 
from the lower order coefficients) have to be added to the right hand side. As the determinant 
of A is non-vanishing except for the three cases mentioned above, we do not have any further 
free parameter than the three parameters from the three linearly independent solutions. 

The solutions with lowest order no = 1 is the one ()4.13p displayed in the main text. The 
most general solution is 



ai{T) 



«2(T) 



768a; sinh^(Ta; 
+(A_3 + A_i + Ai) sinh(4ra; 
1 



4(3A_3 - A_i + 3Xi)Tuj + 8(A_3 - Ai) sinh(2ra;) + 



192a; sinh^(r6< 



4(3A_3 - A_i + 3Ai)ra; cosh(ra;) + 



+(11A_3 + A_i - 9Ai) sinh(ra;) + (3A_3 + A_i - Ai) sinh(3ra;) 



(B.17) 



(B.18) 



192a; sinh^(Ta;) 

+3(5A_3 + A_i - 3Ai) sinh(2ra; 



4(3A_3 - A_i + 3Ai)Ta; + 2(3A_3 - A_i + 3Ai)Ta; cosh(2ra;) + 

(B.19) 



These coefficients describe Hamilton's principal function corresponding to a perturbation term 
(as this is still first order in A) 



A / A_ 



4! 



( ^x\t) + ^x\t)x\t) + \ix\t)] 
\ o;^ a;^ / 



(B.20) 

-3, 



added to the Lagrangian for the harmonic oscillator. The solutions that start with no 
no = — 1 correspond to terms and x^x"^ added to the Lagrangian. 

Applying the same techniques to the last recursion equation ()4.17p reveals that there is only 
one free parameter ^ for the solutions, which does appear in the solution (|4.18p . 



To summarize, we have found all the solutions to the fixed point equations under the (phys- 
ically justifiable) assumption that all the solutions can be expanded into a power series in T 
starting with some lowest, not necessarily positive, order. The free parameters that do appear 
in the solutions can be all interpreted in terms of 'large scale' physics. Hence there are no 
discretization ambiguities left, if we require reparametrization invariance to hold. 
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C Recursion relations for the Pi 

Here we want to shortly explain how one can tackle the recursion relations (|4.14p 

h / 3tanh(Tuj) (n),r^s 3tanh(Ta;) („),^, tanhfTo;) (n) ,r^~. 

M"«'(2T, = ^<"' (T, -J^ + M"' (T) ^ (C.2) 

;i/3tanh(ra^ („) 3tanh(ra;) 
+ ^(cosh2(r.) "° (^)+2cosh(ru;)«i (^^J 

These recursion relations have a homogeneous part 

Pr >i2T) = r (T) (l + ^-1^) + M"' m (C,3) 

,r"(2T) = (D -J^ + (T) -J^ . (C.4) 

to which an inhomogeneity is added. Hence there will be at least one ambiguity, as to every 
solution of the inhomogeneous equations one can add an arbitrary multiple of the solutions to 
the homogeneous solutions. One therefore needs to determine the solutions to the homogeneous 
equations, if one wants to find the full space of solutions. These solutions can be obtained 
as explained in the main text, or alternatively, by subsequently transforming the variables to 
simplify and decouple the equations. Another general method is the expansion in a power series 
in T, which replaces the 'coupling functions' (3o{T), (3i{T) by infinitely many coupling constants 

Assume we are being given the solution /3q , (3^ to the homogeneous equations, in our case 

Pl^{T) = ncoth{Tuj) - V- ^ 



sinh^(ra;) 



/j, , 2TaJcosh(T) 



We want to introduce new variables /3o, /3i such that the recursion relations (|C.ip decouple 
and simplify. To this end we can use the solutions (|C.5p . To decouple the equations we choose 
the new variables such that /3o ~ ^ and /3i ~ for the solutions in the new variables: 

/3o(T) := 2cosh(Tw)/3o(T) + /3i(T) 

/3i(T) := /3o(T) + ^ cosh(ra;)/5i(r) . (C.6) 

Indeed we find for the (homogeneous) fixed point conditions 

/3o*(2T) = 2cosh(Tw)/3o*(T) 

^1{2T) = 2/3^(T) . (C.7) 
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The homogeneous solutions in the new variables are now 

/3^(T) = 2;Usinh(Tw) (C.8) 
/3f (T) = uTu: , (C.9) 

To write the fixed point equations into the form f(2T) — f(T) = g{T) with g{T) representing 
the inhomogeneous terms, we apply another transformation: 

Po{T) := (CIO) 
smh(Ta;j 

^i(r):=M^ . (C.ll) 

This will lead to the following inhomogeneous fixed point conditions (where we replaced the a^""* 
by the fixed point KUh ) 



fiQ{2T) - /3*{T) = ^^ ^ . ^^^^^ ^ cosh(ra;)(10ra; + 2Tuj cosh{4:Tuj) - 4sinh(2ra; 



/3i (2r) - PI (T) = ^ . ^3^^^ ^ - 8Tuj - QTu cosh(rw) - 2Tuj cosh(2Tw) + 



h 

32a;2sinh*(2Ttj) 

sinh(4Ttj)) + 2(-6Ttj cosh(2Tw) + sinh(2ra;) + sinh(4Ttj)) 

32ra;2sinh3(2ra;) 

2rtj cosh(3Ta;) + sinh(TLij) + 3sinh(2ra;) + sinh(3Ta;) + sinh(4rcj 

(C.12) 

The fixed point condition for both recursion relations is now of the form /(2T) — f{T) = g{T) 
where g{T) is some known function of T. For this one can use an 'integration table', see below. 
Alternatively a power series expansion of the solution can be easily given: li g = Y2n^o 9^"^"^ 
then 

/(^) = E^^^" • (C.13) 

Note that it would be inconsistent to have a constant term ~ T° in g{T). An arbitrary integration 
constant can be added to the solution f{T). This ambiguity in the solutions corresponds however 
to adding an arbitrary multiple of the homogeneous solutions (which here are the constants) to 
one solution of the inhomogeneous recursion relations. 
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The following table can be used to obtain a more explicit solution for f{T). 





g{T) = f{2T) - f{T) 


T 


T 


coth^ T 


cosh 21 ' 1 
sinh^ 2T 2 sinh^ T 


cothT 


i 

sinh 2T 


T 

sinh^ T 


OT^ coth 21 
sinh2T 


tantiT 


tanh 1 
cosh 2T 


TcothT 


T tanh T 


T 

sinh 2T 


OT^ 1— cosh 1 
sinh2r 


T tanh T 


rtanh2r + r^ 


In sinh T 


ln(2 cosh T) 


— ^ cothT 




TfiT) 


Tg'{T) 



(C.14) 



In this way one can find all solutions to the fixed point equations for the /3i(T): 



0^(T) = ucothfTw) — uTuj K + f 2 + cosh^fTo;) — 3TtJ cothfTtj' 

°^ ^ ^ ^ sinh2(rw) 32w2sinh2(rw) ^ ^ ^ ^ ' 



_ 2ra;cosh(ra;) _ 2 
^^^^""^ sinh2(rc^) ^'cu^sinMTu;) + 



32sinh^(ra;) 



4Tuj + 2ra; cosh(2ra;) - 3 sinh(2Ta;)) . (C.15) 
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